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Abstract. The precise analysis of properties of planetary nebulae is strongly dependent on good models for the stellar ionizing 
spectrum. Observations in the UV - X-ray wavelength range as well as NLTE model atmosphere calculations of spectra of 
their exciting stars have shown that neither blackbody fluxes nor "standard" NLTE atmosphere models which are composed out 
of hydrogen and helium only are good approximations. Strong differences between synthetic spectra from these compared to 
observed spectra at energies higher than 54 eV (He II ground state) can be ascribed to the neglect of metal-line blanketing. 
Realistic modeling of the emergent fluxes of hot stars in the UV - X-ray wavelength range requires metal-line blanketed NLTE 
model atmospheres which include all elements from hydrogen up to the iron-group. For this purpose, we present a grid (solar 
and halo abundance ratios) of metal-line blanketed NLTE model atmosphere fluxes which covers the parameter range of central 
stars of planetary nebulae. 
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1. Introduction 

Post-AGB stars display the hottest stage of stellar evolution just 
before they enter the white dwarf cooling sequence. Effective 
temperatures up to 180kK have been determined (e.g. Werner 
et al. 1997) by means of spectral analyses based on NLTE (non 
local thermodynamic equilibrium) model atmosphere tech- 
niques. 

A sufficient condition for the use of LTE models in spec- 
tral analyses is, that the collisional rates between atomic states 
are much larger than the corresponding radiative rates (cf. 
Kudritzki & Hummer 1990). This is e.g. fulfilled in the case 
of "cool" (r eff ;40kK) DA white dwarfs (Napiwotzki 1997). 
However, there are always NLTE effects in any star which are 
easily detectable if high-resolution high-S/N spectra, especially 
of the high energy wavelength range (UV, EUV, X-ray), are 
evaluated (e.g. Rauch 1997 and references therein) and these 
effects are particularly important for hot stars like e.g. central 
stars of planetary nebulae (CSPN). 

Observations in the UV - X-ray wavelength range have 
shown that this region is generally dominated by metal opac- 
ities. Since the flux maxima of hot stars lie in the EUV/X- 
ray wavelength (Fig. 1), and the flux in this region is crucial 
to the ionization of ambient matter, e.g. a surrounding plane- 
tary nebulae (PN), the employment of realistic fluxes calculated 
from fully line-blanketed NLTE model atmospheres as ionizing 



Send offprint requests to: T. Rauch 

Correspondence to: Thomas.Rauch@sternwarte.uni-erlangen.de 



spectra in photoionization models is very important. However, 
blackbody fluxes are still often in use. 

In the early 80's of the last century, the calculation of 
NLTE models was rudimental and hampered by the existing 
insufficient numerical methods and computational capacity (cf. 
Rauch 1997). Following Pottasch et al. (1978), who has shown 
that, based on intermediate band photometry within 1500 - 
3300 A, some CSPN have their radiation well represented by 
blackbodies (and others not so well), Kaler (1983) addressed 
that the major problem in PN analyses is the model selection. 
He invoked then that a blackbody has at least the virtue of sim- 
plicity. 

The implementation of the "accelerated lambda iteration" 
(ALI) technique in the NLTE model-atmosphere codes (Werner 
1986) provided an efficient method to calculate realistic syn- 
thetic spectra of hot stars. However, at that time it had been 
impossible to calculate extensive grids of fully line-blanketed 
NLTE models due to the lack of reliable atomic data (cf. Rauch 
2003). Thus, Kwok (1994) repeated Kaler's (1983) statement, 
that due to the variety of CSPN parameters, the blackbody as- 
sumption is probably not worse than any other specific set of 
models. In order to demonstrate the strong deviations at en- 
ergies higher than 13.6eV and 54.4 eV (H i and He n ground 
state thresholds, respectively), we show in Fig. 1 a comparison 
of blackbody fluxes to NLTE model atmosphere fluxes. The 
weaknesses of a blackbody are obvious. Even in a H+He model 
with r e ff = 1 lOkK, log g — 7, and a solar He/H ratio, the model 
flux at energies higher and lower than the He n ground state 
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Fig. 1. Top: Comparison of NLTE model atmosphere fluxes (T e ff = HOkK, logg = 7) calculated with different chemical composition. Short 
dashes: pure H model, thick line: H+He, thin line: H-Ni with [Z] = 0. Note that the main deviations of the blackbody (long dashes, 110 kK) to 
the synthetic spectra occur shortwards of the indicated H I and He II ground state thresholds. Lower panel: Detail of the upper panel; in addition, 
a H-Ni model with [Z] = -1 (thin line), and three blackbody fluxes (long dashes, 90, 110, 130 kK) are shown. The prominent He II and O v 
absorption edges are indicated. Note the drastic decrease of the flux level at shorter wavelengths, if the opacities of metals are considered. 



cannot be matched simultaneously with one blackbody flux. 
Instead, blackbody fluxes with 90 kK and 130kK are needed, 
respectively. At higher energies {A < 160 A), the situation gets 
even worse. 

NLTE model atmosphere fluxes have been successfully 
used, e.g. by Rauch et al. (1994, 1996) for the construction of 
consistent models of PN and CSPN in the cases of the PNe 
K 1-27 and LoTr 4. This approach — consistent models of PN 
and CSPN — appears to be a promising way to improve PN 
analyses. 

Recently, Armsdorfer et al. (2002) selected the spatially 
resolved, round PN NGC 2438 as a test example to check 
for effects of CSPN models on PN shell modeling. For 
their study they used the grid parameters r ei f = 110kK and 



logg = 7 which are closest to the result of Rauch et al. (1999, 
Teff =114+ lOkK, log g- 6.62 + 0.22) from spectral analysis 
of the CSPN of NGC 2438. They found that the use of a black- 
body as ionizing source leads e.g. to an underestimation of the 
[O m] A 5007 A line and an overestimation of the He n A 4686A 
line strength and, hence, an overestimation of the CSPN tem- 
perature, compared to the use of a model from the H-Ca grid 
(Sect. 1). 

Rauch (1997) has investigated in detail on the impact of 
light metals (Li - Ca) on NLTE model atmospheres of com- 
pact hot stars (see this paper for a more detailed descrip- 
tion of advances in NLTE calculations and the consideration 
of metal opacities). He found that CNO opacities reduce the 
flux strongly at energies higher than the oxygen absorption 
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Fig. 2. Comparison of NLTE model atmosphere fluxes (r e ff = 130kK, 
logg = 7, solar abundances) without (thick line) and with consider- 
ation of iron-group opacities. For clarity, the spectra are convolved 
with a Gaussian of FWHM = 0.5 A (the original H-Ni flux is shown 
in addition, shifted down by 0.5 dex). 
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Fig. 3. Parameters of available NLTE model atmosphere fluxes. The 
H-Ni models are limited to r cff < 190 kK because Kurucz's line lists 
(1996) provide only data up to ionization stage ix. 



edges (Fig. 1) while the opacities of light metal reduce the 
flux even more, but at higher energies. Rauch concluded that 
the metal opacities have a great influence on the structure of 
model atmospheres. For the calculation of realistic emergent 
fluxes in the EUV and X-ray range, the consideration of met- 
als is indispensable. Consequently, a first grid of models which 
are composed out of all elements from hydrogen to calcium 
had been calculated with solar ([X] 1 = [Y] = [Z] = 0) and 
halo ([X] = [Y] = 0, [Z] = -1) abundance ratios. This al- 
lows to interpolate or extrapolate slightly in the metal abun- 
dances [Z]. The emergent fluxes are available in tabular form 
at http : //astro . uni - tuebingen . de/~rauch. 



1 [ ]: log (abundance / solar abundance) 



Table 1. Statistics of model atoms used in the calculation of the 
T e ff = 1 lOkK, log g = 7 models. NLTE is the number of levels treated in 
NLTE, RBB (radiative bound-bound) is the number of line transitions. 
495 additional levels are treated in LTE. 



ion 


NLTE 


RBB 


ion 


NLTE 


RBB 


H i 


15 


105 


Mg vi 


4 





H ii 


1 


_ 


Mg vn 


1 





He i 


5 


3 


Ahv 


1 





Hen 


16 


84 


Alv 


5 





He in 


1 


_ 


Alvi 


5 





C m 


3 


1 


Al vn 


1 





C iv 


57 


295 


Sirv 


1 





C v 


1 





Si v 


4 





N iv 


3 





Si vi 


2 





N v 


45 


185 


Si vn 


1 





N vi 


1 





P V 


3 





O iv 


4 


2 


P VI 


1 





O v 


5 


1 


P VII 


2 





O vi 


1 





P VIII 


1 





F iv 


2 





S vi 


3 





F v 


3 





S vn 


6 





F vi 


3 





S vm 


1 





F vn 


2 





CI vi 


3 





F vm 


1 





CI vn 


3 





Ne iv 


4 





CI vm 


1 





Ne v 


9 


3 


CI ix 


1 





Ne vi 


6 


5 


Ar vi 


4 





Ne vn 


4 


2 


Ar vn 


3 





Ne vm 


1 





Ar vm 


3 





Na iv 


1 





Arix 


1 





Na v 


4 





K vi 


1 





Na vi 


5 





K vn 


3 





Na vn 


4 





K vm 


4 





Na vm 


1 





K ix 


3 





Mg iv 


1 





K x 


1 





Mg v 


5 





total 


286 


676 



Table 1. Continued. Numbers in brackets denote individual levels and 
lines used in the statistical NLTE line-blanketing approach for the el- 
ements Ca - Ni. 
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The influence of the iron-group elements was demonstrated 
e.g. by Dreizler & Werner (1993), Lanz & Hubeny (1995), and 
Deetjen et al. (1999). Thus, parallel to the calculation of the H- 
Ca grid, test calculations have been performed in order to check 
the impact of the iron-group opacities (Rauch 2000, Rauch & 
Deetjen 2001, Rauch 2002). The consideration has an influence 
on the atmospheric structure and thus, also on the strengths of 
the absorption edges of light metals. Since significant changes 
in the flux level in the UV - X-ray region occur (Fig. 2), the 
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Table 2. Example header of a flux table. 

* TUEBINGEN NLTE MODEL ATMOSPHERE PACKAGE 

* 

* DATA FILE CREATED BY 02W AT 19.11.02 11:47:43 

* MODEL CALCULATED BY PR02 AT 20.01.01 03:59:57 

* 

* T eff = 110000 K 

* log g = 7.00 (cgs) 

* 

* T eff from integrated flux: 110599 (quality check) 

* 

* NORMALIZED NUMBER FRACTIONS 

* 

* H 9.079E-01 

* HE 9.O79E-02 

* C 3.614E-04 

* N 9.079E-05 

* 7.211E-04 

* F 3.614E-08 

* NE 3.614E-05 

* NA 1.811E-06 

* MG 2.871E-05 

* AL 2.281E-06 

* SI 3.614E-05 

* P 2.281E-07 

* S 1.439E-05 

* CL 1.811E-07 

* AR 4.55OE-06 

* K 1.143E-07 

* CA 4.086E-05 (GENERIC MODEL ATOM: CA - NI) 

* 

* COLUMN 1: WAVELENGTH GIVEN IN A 

* COLUMN 2: ASTROPHYSICAL FLUX ( F = 4H ) 

* F LAMBDA GIVEN IN ERG/CM**2/SEC/CM 



* 19951 LINES FOLLOWING 
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Fig. 4. Comparison of NLTE model atmosphere fluxes (log g - 7) cal- 
culated with solar (top) and halo (bottom) abundance ratios at three 
different temperatures. 



iron-group opacities should not be neglected. Thus, a second 
grid with all elements from hydrogen up to the iron group has 
been calculated. This is described in the following. 

2. NLTE model atmospheres 

Our NLTE models (Werner 1986, 1989; Werner & Dreizler 
1999, Werner et al. 2003) are plane-parallel and in hydrostatic 
and radiative equilibrium. Hence, our models neglect wind ef- 
fect onto the spectrum, which can be significant close to the 
Eddington limit (Gabler et al. 1989). 

All elements from hydrogen to the iron-group elements 
can be considered simultaneously. More than 300 levels (lim- 
ited by numerical accuracy in the used code-version) can be 
treated in NLTE with more than 1 000 line transitions cal- 
culated in detail. Millions of lines of the iron-group ele- 
ments tabulated in Kurucz (1996) and data from the Opacity 
Project (Seaton et al. 1994) are accounted for using a statisti- 
cal approach and an opacity sampling method (Deetjen 1999, 
http : //astro. uni-tuebingen. de/~deetj en). 



For the construction of the model atoms used in our cal- 
culations, we follow Rauch (1997). H, He, C, N, O, and Ne 
are represented by small model atoms, the other species are 
modeled to fit into the capacity limit of the code. As a conse- 
quence, many of the ions of F, Na - K, (Table 1) are too small to 
include radiative transitions, not even the resonance lines. We 
performed some test calculations with larger model ions for se- 
lected atoms in order to include at least all resonance lines of 
these ions. We found that these lines are very weak, so they do 
not have any significant influence on the atmospheric structure. 
Since our aim is not a precise spectral analysis of a given star 
but the calculation of emergent fluxes as ionizing spectra, this 
is a good approximation. 

In contrast to Rauch (1997), Ca is not treated as a sin- 
gle species but combined with the iron-group elements into a 
generic model atom in order to consider all Ca lines given by 
Kurucz (1996). Within the generic model atom, Ca - Ni are 
represented with solar abundance ratios. Depending of r ei f , the 
model atoms have individually been adjusted due to the dom- 
inant ionization stages. The statistics of representative model 
atoms are summarized in Table 1 . For a more detailed descrip- 
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tion of the used atomic data and the consideration of the iron- 
group elements, see Rauch & Deetjen (2003). 

We calculated a grid of models within r e fj =50-1 OOOkK 
and logg = 5-9 (Fig. 3), with solar and halo abundances 

(Fig. 4). We use r™ d = (4n / cr £° <H v dv) 1/4 for the conver- 
gence criterion (^H v is the Eddington flux). If this is closer than 
2% to the target r e ff, the model calculation is stopped in order 
to stay within a reasonable cpu time on the CRAY computers. 
We mention that it may take about one day to calculate a single 
model. 

Flux tables have subsequently been calculated within 
5-2000A binned to 0.1 A intervals. r™ od as well as the 
model's T e s and logg are given in the tables' head- 
ers (Table 2). These tables are directly suited to be used 
with the photoionization code CLOUDY (Ferland 2000, 
http://www.pa.uky.edu/~gary/cloudy/, versions C96 
Beta 4 and later). The complete flux tables, based on the fre- 
quency grids used in the model calculations, are also available 
in case that a normalization to a certain (e.g. V) magnitude is 
required. 

In order to judge the necessity of a detailed consideration 
of iron-group opacities in calculations of ionizing fluxes which 
shall be used in photoionization models of a PN, we measured 
the changes of the strengths of the H i and He n absorption 
edges in the H-Ca and H-Ni fluxes. We calculated 



s = 



Aq-AA A 



Ja -aa a 



c *dA 



r-lo 



(1) 



with Ao = 91 1.5 A and Aq = 227.8 A for H i and He n, respec- 
tively, and AA = 50 A. While the changes of the H i absorption 
edge are generally less than 10%, the strengths of the He n ab- 
sorption edge may change by a factor of about 2 - 3 at lower 
temperatures. Depending on \ogg, there is a maximum of s, 
e.g. s = 2.3 at T e g = 70kK and log g = 5. s decreases towards 
lower and higher T e s and is about 1 ±0.2 at T e s > 120 kK. 

A similar grid for hydrogen-deficient stars (Fig. 5) has 
been calculated. For this grid we selected the "typical" abun- 
dance pattern (He:C:N:O=33:50:2:15 by mass) of PG1159 
stars (Werner 2001). This allows to estimate abundance effects 
onto the emergent flux. 

3. Conclusions 

Presently available NLTE model atmosphere codes have ar- 
rived at a very high level of sophistication (Mihalas 2003, 
Werner et al. 2003). The accelerated lambda iteration (ALI) 
method represents a powerful tool to calculate realistic metal- 
line blanketed atmospheres in a reasonable time, even on nowa- 
days available PCs — these got such a tremendous number- 
crunching power in the last two years that they can now easily 
compete with CRAY and other (super)computers. The statisti- 
cal approach which has been used to construct a generic model 
atom for all elements from Ca to Ni (Sect. 2) has been improved 
and can now be employed to construct model atoms for other 
light metals in order to consider all of their lines, provided by 
Kurucz (1996) and the Opacity Project (Seaton et al. 1994), in 



detail. Thus, there is no need to reflect on blackbodies or sim- 
ple model atmospheres as ionizing sources for photoionization 
models. 

To summarize, the use of metal-line blanketed NLTE model 
atmosphere fluxes, which consider at least the CNO opacities in 
detail, as ionizing spectra for photoionization models is highly 
recommended. The here presented fully metal-line blanketed 
NLTE model atmosphere fluxes which consider all elements 
from hydrogen up to the iron group are well suited for this pur- 
pose. 

However, the user of these model atmosphere fluxes should 
be aware that they should only be used in a PN analysis if 
the CSPN abundances are in agreement with the model abun- 
dances. The use of fluxes from models with solar abundance 
ratios will spoil any PN analysis if the CSPN is hydrogen de- 
ficient because the flux shortwards of the He n ground state 
edge can be different by about one dex (Fig. 5). A preceding 
spectral analysis of the CSPN which yields r e ff, logg, and 
abundance ratios and an individually calculated small "grid" 
of NLTE model atmosphere fluxes may improve a subsequent 
PN analysis — and allows to construct consistent models of 
PN and CSPN (Sect. 1). This is also a prerequisite for reliable 
hydrodynamical models of the PN. 
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